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We study theoretically and experimentally the response of a microwave superconducting stripline 
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A theoretical model which attributes the self modulation to a thermal instability yields a good 
agreement with the experimental results. 
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I. INTRODUCTION 

Nonlinear effects in superconductors have significant 
implications for both basic science and technology. 
Strong nonlincarity may be exploited to study some im- 
portant quantum phenomena in the microwave region, 
such as quantum squeezingJ^i^ and experimental obser- 
vation of the so called dynamical Casimir effect^. These 
effects may also allow some intriguing technological ap- 
plications such as bifurcation amplifiers for quantum- 
limited measurements^!^, resonant readout of qubitsi, 
mixers^, single photon detectors^, and more. 

In this paper we study theoretically and experimen- 
tally the response of a superconducting (SC) microwave 
stripline resonator, designed for enhanced nonlincarity, to 
a monochromatic injected signal. We find that there is a 
certain range of driving parameters, in which a novel non- 
linear phenomenon immerges, and self-sustained modu- 
lation (SM), of the reflected signal from the resonator, 
is generated. That is, the resonator undergoes lim- 
ited cycle oscillations, ranging between several to tens 
of megahertz. Similar phenomenon was briefly reported 
in the 60'9i2iiiii^ii^ in dielectric resonators which were 
partially coated by a SC film, but it was not thor- 
oughly investigated and its significance was somewhat 
overlooked. This phenomenon is of a significant impor- 
tance as it introduces an extreme nonlincarity, which is 
by far stronger than any other nonlincarity observed be- 
fore in SC resonatorsii. It results in a very strong in- 
termodulation gain, strong noise squeezing and period 
doubling of various orders^^, strong coupling between dif- 
ferent resonance modest, and more. 

The central results of this study have been recently re- 
ported in a short paperiS,. In this paper we extend the 
previous report and derive a theoretical model, accord- 
ing to which the SM originates by a thermal instabil- 
ity in the resonator. The numerical integration of the 
model's equations of motion exhibits SM, which have 
similar characteristics as the experimental results. In 
addition, we derive analytic expressions for the expected 
SM frequency and the spectral power density. These ex- 
pressions also yield a good agreement with the experi- 



mental results. 

This paper is organized as follows. First we briefiy 
describe the design of our devices and the experimental 
setup. Then we present the SM phenomenon, as mea- 
sured in our devices. Afterwards, we derive the theoreti- 
cal model, discuss and justify its underlying assumptions 
and present the numerical integration results. Finally, 
we quantitatively compare the predictions of the model 
to typical experimental results. 

II. EXPERIMENTAL SETUP AND CIRCUIT 
DESIGN 

The majority of the experiments are performed using 
the experimental setup described in Fig. [Tfa). The res- 
onator is stimulated with a monochromatic pump tone at 
an angular frequency Wp. The power refiected off the res- 
onator is amplified at room temperature and measured 
by using both a spectrum analyzer (SA) in the frequency 
domain, and an oscilloscope, tracking the reflected power 
envelope in the time domain. In other experiments, the 
jS'iil reflection coefficient is measured using a network 
analyzer (NA), connected directly to the resonator RF 
port. All measurements are carried out while the device 
is fully immersed in liquid Helium. 

A simplified circuit layout of the device is illustrated 
in Fig. [Ijb). The resonator is designed as a stripline 
ringl^iii, having a characteristic impedance of 50 Q,. It 
is composed of Niobium Nitride (NbN) deposited on a 
Sapphire wafer. The first few resonance frequencies fall 
within the range of 2 — 8 GHz. A feedline, weakly coupled 
to the resonator, is employed for delivering the input and 
output signals. A microbridge, which is employed as a 
weak-link that allows the manipulation of the resonator's 
resonance frequencies^^, is monolithically integrated into 
the structure of the ring. Its angular location, relative 
to the feedline coupling location, maximizes the RF cur- 
rent amplitude flowing through it in one of the resonance 
modes, and thus maximizes its coupling to that mode. 
Further design considerations, fabrication details as well 
as normal modes calculation can be found elsewhere^^. 

The results presented herein are obtained using two 



2 




Oscilloscope 



Spectrum Analyzer 



NbN 



FIG. 1: (a) SM measurement setup, (b) Schematic layout of the device. Subplot (c) and (d) exhibit optical microscope images 
of the straight and meander shaped microbridges respectively. 
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FIG. 2: (Color Online) |Sii| reflection measurements taken with (a) E15 and (b) E16 devices. Note that the |Sii| reflection 
coefficient is only defined for the case of a steady state refiection from a device. Therefore the 15*111 measurements taken at the 
SM region should be interpreted as an average over time of the l^iil coefficient. 



distinct devices, labeled as E15 and E16, which differ 
by the geometry of the microbridges and by their thick- 
nesses. El 5 has a 1 X 10 //m^ microbridge geometry 
(Fig. lUc)) and a thickness of 200 nm, whereas E16 has 
a 4 X 4 /im^ meander shaped microbridge geometry (Fig. 
[TJd)) and a thickness of 8nm. The meander consists of 
nine strips, where each strip has a characteristic area of 
0.15 X 4 /im^ and the strips are separated one from an- 
other by approximately 0.25 /inJ^. 

The difference in the frequency response between E15 
and E16 can be observed in a simple 15*111 reflection mea- 
surement, obtained using a NA. Fig. [5] shows various 
15*11 1 curves as a function of the pump frequency central- 
ized on the third resonance frequency /a = 5.666 GHz 
of E15 (panel (a)) and the second resonance frequency 
/2 = 3.87 GHz of E16 (panel (b)), labeled by na and m2 
respectively. Each curve represents a measurement with 
a different pump input power. For clarity, the curves 
are vertically shifted upwards, for increasing power val- 
ues. The anomaly of the response is described as follows. 



Above some power threshold the 15*111 line-shapes cease 
having a normal Lorenzian shape, their values substan- 
tially increase, and the resonance curves substantially 
broaden and have steep edges. This behavior contin- 
ues and intensifies as the pump power further increases, 
until eventually no resonance is detected. Furthermore, 
as seen in panel (b), at relatively high power levels, the 
resonance curve of E16 is reconstructed at a new reso- 
nance frequency, red shifted by approximately 45 MHz 
relative to f2- The Lorenzian line-shape of the new res- 
onance frequency, labeled by represents a linear be- 
havior in that power range. These experimental results 
suggest that external stimulation can cause a significant 
resonance shift in E16, while E15 can only experience an 
increase in its damping rate. 
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FIG. 3: (Color online) Typical experimental results of the SM phenomenon in the frequency domain. Panel (a) plots a colormap 
of the reflected power Prefl as a function of the pump power Ppump and the measured frequency fsA centralized on the third 
resonance frequency fs {fg\ — fsA — fa) of E15, while the resonator is stimulated with a monochromatic pump at /s. Panel 
(b) plots a colormap of the reflected power as a function of the centralized measured frequency around the pump frequency 
fsA ~ fsA — L0p/2TT and the centralized pump frequency around the resonance frequency /pump ~ i^p/2n — /s. The resonator 
is stimulated by a monochromatic pump having a power of Ppump = —29.35 dBm which drives the resonator into the regular 
SM zone. 



III. SM EXPERIMENTAL OBSERVATION 

We now turn to investigate the region where SM im- 
merges. Fig. [3] shows typical experimental results of the 
SM phenomenon in the frequency domain, as measured 
with E15. The dependence of the SM on the pump power 
is shown in panel (a) and described as follows. At low 
input powers, approximately below —33.25 dBm, and at 
high input powers, approximately above —25.0 dBm, the 
response of the resonator is linear, namely, the reflected 
power from the resonator contains a single spectral com- 
ponent at the frequency of the stimulating pump tone Up. 
In between the two linear regions, there is a rather large 
power range in which regular SM of the reflected power 
from the resonator occurs. It is realized by rather strong 
and sharp sidebands, which extend over several hundreds 
megahertz at both sides of the resonance frequency. The 
SM frequency, which is defined as the frequency difference 
between the pump and the primary sideband, increases 
as the pump power increases. 

The regular SM starts and ends at two power thresh- 
olds, referred to as the first and the second power thresh- 
olds. The first power threshold occurs at a very narrow 
power range of approximately 10 nW, during which the 
resonator's response desists being linear. It experiences 
a strong amplification of the noise floor, known also as 
noise rise, over a rather large frequency band, especially 
around the resonance frequency itself. The second power 
threshold occurs on a slightly larger power range than the 
first one and has similar, but less extreme characteristics. 

As shown in panel (b), the dependence of the SM 
on the pump frequency is rather symmetric around the 



resonance frequency. It occurs only within a well de- 
fined frequency range around the resonance frequency. 
A small change in the pump frequency can abruptly ig- 
nite or quench the SM. Once started though, the mod- 
ulation frequency has a relatively weak dependence on 
the pump frequency. E16 exhibits similar behavior, with 
slightly different properties, as discussed below in section 
IIVA3bl 



IV. THERMAL INSTABILITY 

In this section we propose a theoretical model accord- 
ing to which the SM originates by a thermal instabil- 
ity in the SC stripline resonator. Current-carrying su- 
perconductors are known to have two metastable phases 
sustained by Joule self-heating2S. One phase is the SC 
phase and the other is an electrothermal local phase, 
known as hotspot, which is basically an island of normal- 
conducting (NC) domain, with a temperature above the 
critical one, surrounded by a SC domain. This phe- 
nomenon can be explained by the heat balance equa- 
tion holding at more than one temperature. Due to an 
extcrnaP^ or internaP^ perturbation, the hotspot can re- 
cover to the SC phase or vice versa and thus oscillates be- 
tween these phases. Such self-sustained oscillations were 
often observed in experiments, for the case of a SC mi- 
crobridge driven by an external dc voltage or current (see 
review^S and references therein). 

In the present the microbridge is integrated 

into a stripline resonator, the system is driven into in- 
stability via externally injected microwave pump tone. 
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FIG. 4: (Color online) Schematic model of the driven res- 



GuJoS {t — t'). Consider the case of relatively high tem- 
perature, hgT^g 3> hujQ, where fc^ is Boltzmann's con- 
stant, and Teff is a weighted average between T and Tq, 
where Tq is the temperature of the coolant and the weight 
factors are discussed in Ref."'. Then at steady state, the 
variance of the noise is given by G = 2jkBTcg/hu!Q. 

The steady state solution of Eq. ([1]), which is denoted 
as -Boo, is given by 

J (wp - Wo ) - 7 



2. Thermal Balance 



Nonlinearity, according to our simple theoretical model, 
results from the coupling between the equation of mo- 
tion of the mode amplitude in the resonator (Eq. ([T])), 
and the thermal balance equation (Eq. ([3])) in the mi- 
crobridge. The mechanism presented in this model is 
somewhat similar to one of the mechanisms which cause 
self oscillations in an optical parametric oscillator—. 



A. Equations of Motion 

This section presents the equations of motion of the 
mode amplitude in the resonator and the thermal balance 
in the microbridge. Note that Eqs. (P)-®, (l5|)-(I10l) also 
appear in Ref^, and are re-presented here to ensure that 
the paper is self-contained. 



1. Mode Amplitude 

Consider a resonator that is driven by a weakly cou- 
pled feedline carrying an incident coherent tone a'" = 

{\b" 

ujp is the driving angular frequency (See Fig. 2]). The 
mode amplitude inside the resonator can be written as 
A = i?(t)e~*"p', where B (t) is a complex amplitude 
which is assumed to vary slowly on a time scale of 1/wp. 
In this approximation, the equation of motion of B reads^ 



''p , where 5'" is a constant complex amplitude 
oc Ppump, where Ppump is the driving power) and 



dB 



= [z (c.p-c.0 (T))-7(r)]i?- 



(1) 



where coq (T) is the temperature dependant angular res- 
onance frequency, T is the temperature of the hotspot, 
7 (T) = 7i -|- 72 (T), where 71 is the coupling constant 
between the resonator and the feedline, and 72 (T) is the 
temperature dependant damping rate of the mode. 

The term c™ represents an input Gaussian noise 
with a zero-mean and a random phase, thus (c™) — 
0, (c'"(i)c'"(i')> = (c'"*(i)c"*(t')) = 0, and its au- 
tocorrelation function is given by {d'^[t)c!-'^* (t'')) = 



Consider the case where the nonlinearity originates by 
a local hotspot in the resonator's microbridge. If the 
hotspot is assumed to be sufficiently small, its tempera- 
ture T can be considered as homogeneous. The temper- 
ature of other parts of the resonator is assumed be equal 
to that of the coolant Tq. The power Q heating up the 
hotspot is given by Q = ^Qt, where Qt = fiwo272 \Bf is 
the total power dissipated in the resonator, and ^ k ^ 1 
represents the portion of the dissipated power that is be- 
ing absorbed by the microbridge. The heat balance equa- 
tion reads 



AT 

C— = Q (B) - IF, 



(3) 



where C is the thermal heat capacity, W — H {T — Tq) 
is the heat transfer power to the coolant, and H is the 
heat transfer coefficient. 

The steady state solution of Eq. ([3]), which is denoted 



as Too, f^or B 



Boo, is given by 



To 



Kfkjo2^2 \Bo 
H 



(4) 



3. Stability zones 

a. Coupling mechanism The coupling mechanism 
between Eq. ^ and Eq. ^ is based on the dependence 
of both the resonance frequency and the damping rate of 
the driven mode on the resistance of the microbridge^?^, 
which in turn depends on the temperature T— . 

Here we assume the simplest case, where the resonance 
frequency luq, and the damping rate 72 have a step func- 
tion dependence on the temperature of the hotspot (the 
step occurs at the critical temperature Tc). It is based 
on the fact that recent experiments with photodetectors, 
based on a thin layer of NbN, have demonstrated an in- 
trinsic switching time on the order of 30 ps (see Ref.— and 
references therein). In addition, when illuminating our 
devices with a power modulated infrared light in a simi- 
lar way to the experiment described in Ref.^^, we measure 
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FIG. 5: Stability zones of the system, corresponding to (a) 
E15 and (b) E16 related assumptions. The solid lines are 
obtained from Eq. © using ]-Boo|^ ~ Es and |-Boo[^ ~ En- 



diagram. In the monostable zone either the SC phase or 
the NC phase is locally stable, whereas in the bistable 
zones both phases are locally stable^ii^. In the null- 
stable zone, on the other hand, none of the phases are 
locally stable, and the resonator is expected to oscillate 
between these two phases. As the two phases significantly 
differ in their reflection coefhcients, the oscillations are 
translated into a modulation of the reflected pump tone. 
Note that, the stability diagram indicates the existence 
of both power and frequency hysteresis in the system's 
respons o^'^'^^ . Furthermore, panel (b) in fig. jSj shows 
that the dependence of the SM on the pump frequency 
is asymmetric in the case where a resonance frequency 
shift occurs. 

This asymmetry is indeed observed experimentally 
when measuring the SM frequency as a function of the 
pump frequency and power as shown in Fig. [6] panels 
(a) and (b), presenting data obtained with E15 and E16 
respectively. In both cases SM occurs in the nuUstabil- 
ity zone (compare with Fig. [5]). One clearly notices that 
the SM, as measured with E16, is strongly asymmetric in 
frequency in contrast to the case of E15. The maximum 
measured SM frequency is approximately 41.1 MHz and 
57.6 MHz with E15 and E16, respectively. 



B. Adiabatic Approximation 



a clear response up to modulation frequencies of several 
gigahertz with ElO** and several hundred megahertz with 
E15. Thus the transition through the instability point is 
very fast on the time scale of the SM frequency. We fur- 
ther assume that all other parameters in the model are 
temperature independent. 

Under this assumption, the heat generation on the 
temperature T have step-like dependence-*^. The sys- 
tem thus may have in general up to two locally stable 
steady states, corresponding to the SC and NC phases. 
A SC steadj state exists when Too < Tc, or alternatively 
when l-Bool < E^, where Eg = H (Tc — Tq) /2Kj2sf>^a- 
Similarly, a NC steady state exists when Too > 7c, 
or alternatively when |i3oo| > E^, where E^ — 
H (Tc — To) /2K72nfiwo, where the subscript (s) and (n) 
denote the value of the corresponding parameter when 
the system is in the SC and the NC phases respectively. 

b. Stability The stability of each of these phases de- 
pends on both the power and frequency parameters of 
the injected pump tone, as described by the stability di- 
agrams in Fig. |SJ Panel (a) shows a case similar to E15, 
where the SC and the NC phases differ by the damping 
rate value, which has a larger value in the latter phase. 
Panel (b), on the other hand, shows a case similar to E16, 
where the resonance frequency in the SC phase is higher 
than the one in the NC phase, whereas the damping rates 
are the same in both phases. 

Four different stability zones can be identified in the 



In the following section we derive two analytic expres- 
sions, one for the SM period (Eq. (HH)) and another for 
the SM spectral density (Eq. (US]))- Both expressions 
are valid for input powers that are slightly larger than 
the first power threshold. The derivation assumes that 
the system is in the adiabatic regime, namely that the 
rate in which the temperature of the hotspot changes is 
much faster than the SM frequency. In addition it as- 
sumes for simplicity that the pump frequency equals the 
resonance frequency. 



1. Dimensionless variables 

In terms of the dimensionless time r = idot Eq. ([T]) 
reads 



db , , d"" 
dr ajQ 



(5) 



where b = B — Boo, and A = [7 — i (ujp — loq)] /loq. 

Defining the dimensionless temperature O = 
(T — Tq) / (Tc — Tq), and using the dimensionless 
time T Eq. ^ reads 



^+5(6-600) = 0, 



(6) 



where 



gC'iTc-To) 



^09 
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FIG. 6: (Color online) Measured SM frequency /am as a function of the pump power Ppump and the pump frequency /pumpi 
centralized on (a) the third and (b) the second resonance frequencies of E15 and E16, respectively. 



is the steady state value of the dimcnsionlcss temperature 
for a fixed mode amplitude, and 



P 



H 

hujQ 



C (Te - To 



(8) 



(9) 



The steady states solution of Eq. ([5]), where B — Boo, is 
denoted as 



gC{T,-To) 



^05 



2. Adiahatic Solution 



where ^ (r) 6 (t) - Q^o (r). Using Eq. (HH) the so- 
lution of Eq. p3|) for the case gr ^ 1 can be written 
as 



n 1 d9oo (r) 



(14) 



Thus the lagging of the temperature Q (r) behind the 
asymptotic value Qao depends on the rate of change of 
\B\^ (see Eq. (0)). 

Moreover, when the pump frequency equals the reso- 
nance frequency, namely cjp = o^o , B is purely imaginary 
and its time evolution is given by 



B (t) - Bo 



exp 



7r 

Wo 



(15) 



where i?o 



B (0) - Bo 

— iv^27i6'"/7. Consequently Eq. (fT4|) reads 



Assuming the case where Q =/= I in the time interval 
(0, r), and disregarding noise, the solution of Eq. ([5]) is 
given by 



B (r) = Bo 



B (0) - Boo , , . 
1 + exp (--At) 

-t>oo 



(11) 



thus, using the notation (3 = {B (0) — Boo) /Boo and Eq. 
(fTO| . one has 



Ooo (r') = eooo[l + /3exp (-At') + 
+ /?* exp (-AV) + 

+ |/3|%xp[- (A + A*)t']. 



(12) 



In the adiabatic limit, where j/gujo ^ 1; one expects 
that the temperature closely follows the evolution of the 
mode amplitude, namely 8 (t) ~ 8oo (t), thus it is con- 
venient to rewrite Eq. ^ as 



dT 



5^ 



d9oc 
dT 



(13) 



e (t) = Boo (r) - [Oooo - Ooo (r)] . (16) 
5^0 



3. Null-stability zone 

Consider the case of operating in the nuUstability 
zone. Switching between SC and NC phases occurs when 
Q (t) — 1. At that time the mode amplitude B can be 
found from the value of Qoo- Using Eq. (|16p one finds 
to first order in ^/gujQ 



1 + 



7(eooo-i) 



gujo 



(17) 



Using Eq. ([7]) and Eq. PH)) one finds that the mode 
amplitude at switching, to first order in 7/(7070, is given 

by 



\B\ = \Bn 



2gujo 



Boo 



(18) 



7 



where |i3o| = '^o5/2k72P is the value of \B\ for which 
&00 = 1- We denote by (B^) the value of |i?| when 
a switching from the SC to the NC (NC to SC) phase 
occurs. 

a. SM Period The SM period is denoted as 
T = Ts+Tn, where % {T^J is the time, in which the 
system is in the SC (NC) phase. Using Eq. (fT5|) one 
finds 



^ 1 1 -^11 Bqqq 

% = — log^ — 

7s Bs - Boos 



— log^ 

7n Bn 



Bo 



Bo 



(19) 



Slightly above the first SM power threshold one has Bg ~ 
Boos- In this case one expects that % » 7^. Moreover 
writing % as 



1 

Is 



log- 



Bo 



B„ 



B, 



log 



Os 



Os 



Boos — B^ 



(20) 



where Bqs — —W'^lisb™/^a, ^o" the input amplitude 
associated with the first SM power threshold, and ne- 
glecting the first term, which is much smaller than the 
second one, and using Eq. yield 



TABLE I: Model's Parameters 





E15 


E16 




E15 


E16 


loo/2ti [GHz] 


5.7 


3.8 


Q-factor 


880 


250 


d [nm] 


200 


8 


9 [xlO-3] 


9.69 


250 


-Ppump [dBm] 


-25.5 


-49 


P [xlO-«] 


0.80 


1180 


l^nlth [dB] 


-4.15 


-13 




197 


27 



where Wfe = ik-K /T. Thus, using Eq. the spectral 

density slightly above the first SM power threshold is 
given by 



271 (gs - Bn) 



(26) 



r 



Ts 

1 

7s 



log ■ 



Bo 



Bqs 



2gujo 



Boo 

'b^ 



(21) 



Thus, using the notation i) = (6'" — fog") /b™ one finds 
that slightly above the first threshold, when d ^ 1, 
namely when regular SM with a relatively long period 
occurs, the SM period is given by 



T ~ — log —, 

7s ^ ( 1 _ ^ 



1 , 1 
— log^. 

7s w 



(22) 



Note that disregarding noise can not be justified very 
close to the first power threshold, since in that region 
the system is extremely sensitive to fluctuations. 

b. Spectral Density The output signal reflected from 



the resonator is written as a° 



''p', where 6°"* is 



a complex amplitude. According to the input-output re- 
lation, which relates the output signal to the input one^^, 
the following holds 



(23) 



Above the first power threshold the amplitude B (t) is 
periodic, B {t) — B {t + T). If the assumption % ^ 
holds, one finds 




B (t) ~ B„ + (Boos - B„) (1 - e 



(24) 



where the time interval in which the hotspot is in NC 
phase is neglected. The power spectrum of the k*"^ har- 
monic of 6°"* is given by 



Pk 



1 

T 



T 



(25) 



Jf.. Validity of the adiabatic approximation 

We now return to the adiabatic approximation 
j/gujQ <C 1, and examine its validity by estimating the 
value of the parameter g in Eq. ([8]). Consider the case 
where the nonlinearity originates by a hotspot of lateral 
area A^s, forming in the microbridge. The heat capac- 
ity C of the hotspot can be expressed as C = CvAcsd, 
where Cy is the heat capacity per unit volume, and d 
is the thickness of the NbN film. By further assum- 
ing that the generated heat is cooled mainly down the 
substrate rather than along the film.^^, the heat trans- 
fer coefficient reads H = aA^s, where a is the thermal 
surface conductance between the NbN film and the sub- 
strate. According to this notation Eq. ([5]) is expressed as 
g — a/Cydujo- To obtain an estimate for the parameter 
p in Eq. we evaluate the total dissipated power Qt 
in the resonator at the first SM power threshold Eg, by 
assuming that the power dissipated in the microbridge is 
given by nQt = k{1 - jS*!! |th)-Ppump, where Ppump is the 
injected input power, |>S'ii|^jj is the reflection coefficient 
at the first SM power threshold and k ~ 1. On the other 
hand, this equals to the heat flow from the microbridge to 
the substrate W = aA^e (T - Tq). Thus Eq. © can be 
expressed as p = afkoo/CydKQt- The value of a and 
parameters as estimated for NbN on a Sapphire substrate 
at temperature T = 4.2 K^ZjIS, are a ~ 12.5Wcm~2K~^ 
and Cy ~ 2.7 x 10^ Jem" -3K"\ The various measured 
and calculated parameters are summarized in table [H 
which shows that the adiabatic assumption is justified 
for E16, and is marginal for E15. But as the above es- 
timation does not take into account the direct contact 
between the sample and the liquid Helium, it is reason- 
able to assume that the adiabatic assumption is justified 
for E15 as well. 
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FIG. 7: (Color online) Numerical integration of the equations of motion of the model at three distinct pump powers (a) the 
first power threshold (b) some intermediate power between the first and second power thresholds , and (c) the second power 
threshold. Each panel has three subplots. Subplot (i) describes the normalized mode amplitude _Bn (blue), and the normalized 
hotspot temperature Tn (red), as a function of time, where both are normalized by their critical values at which a transition 
from the SC to the NC phase occurs. Subplot (ii) plots the reflected power off the resonator Prefl normalized by the impinging 
pump power Ppump. Subplot (iii) describes the normalized steady state solution of the mode amplitude (Eq.((2])), as a function 
of the normalized frequency, for the case where it is decoupled from the bridge temperature. The solution is again normalized 
by its critical value Es, and the frequency is normalized by the SC resonance frequency. The solid and dashed portions of the 
curves represent solutions which are stable and unstable respectively, according to the stability diagram in Fig. [5jb) . The blue 
and red curves are solutions for the cases where the system is in the SC and the NC phase, respectively. The magenta and 
green dashed lines in subplots (i,iii) show Es and En normalized by Es respectively. 



V. NUMERICAL INTEGRATION 

Following the discussion in section [TTl the numerical 
results presented for E15 are calculated by assuming a 
significant increase in the damping rate and a negligible 
shift in the resonance frequency as the critical temper- 
ature is exceeded, and the results presented for E16 are 
calculated by assuming the opposite case, where only loq 
is temperature dependant. 

Fig. [7] shows numerical integration results of the cou- 
pled Eq. IH) and Eq. ([3]), using the case where a shift in 
the resonance frequency occurs (E16). Results obtained 
by assuming a significant increase in the damping rate 
with no frequency shift are presented in Ref.^^. The res- 
onator is stimulated by an impinging pump tone at the 
SC resonance frequency. Panel (a) shows results that 
are calculated for a pump power in the first SM power 
threshold range. If the system would have been noise- 
less then both the mode amplitude and the temperature 
could reach a steady state. At this steady state the re- 
flection from the resonator is relatively low as the pump 
frequency coincides with the resonance frequency. As 
this steady state is on the edge of instability, the thermal 
noise at a temperature of 4.2 K makes the system un- 
stable and it occasionally falls off the edge and switches 
to the NC phase. When this happens, the dissipation 
slightly increases but more significantly, the resonance 
frequency shifts, and consequently the mode amplitude 
starts decreasing. As a result the heat production, which 
is proportional to the mode amplitude squared, decreases 
and thus when the excess heat is transferred to the sub- 
strate, and the temperature of the bridge decreases be- 
low the critical one, the resonator switches back to the 
SC state and a new buildup cycle of electromagnetic en- 



ergy begins. Accordingly, the reflected power is low for 
most of the time, and these cycles are realized as spo- 
radic but correlated spikes of the reflected power. The 
heat generated in such a spike, raises the probability for 
a sequential spike to occur, and thus induces a positive 
correlation between the spikes (bouncing). The dynamics 
of the relaxation cooldown of these spikes is similar, and 
thus their line-shapes also. The power range in which 
these spikes are triggered, and hence the width of the 
power threshold range, is governed by the noise inten- 
sity 

Panel (b) shows results calculated for some pump 
power in the range of the regular SM. The evolution of 
the system is similar to the one just described, with one 
major difference. When the system is in the SC phase 
the mode amplitude is built toward a nuUstable state and 
thus a steady state is not achieved. As a result, regular 
oscillations occur without the assistance of noise, which 
in general, has a negligible impact in that power range. 

Panel (c) shows results calculated for a pump power 
at the second SM power threshold while assuming a 15 K 
thermal noise and a slightly enhanced damping rate due 
to an average increase in the microbridge temperature. 
The behavior of the system at this threshold resembles 
the first one, but the SC and NC phases exchange roles. 
In this power range the resonator is in the NC, high reflec- 
tive phase for most of the time and noise-induced spikes 
temporarily drive it to the SC, low reflective phase. The 
internal thermal noise at the second threshold is stronger 
than at the first one and consequently this power thresh- 
old range is wider. 

Fig. [S] shows the envelope line-shape of the reflected 
power when a regular SM, having a frequency of approx- 
imately 6 MHz, occurs. The experimental data is seen 
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FIG. 8: (Color Online) SM reflected power line-shapes of ex- 
perimental (blue) and the numerical integration of the model's 
equations of motion (dashed red) results, normalized by the 
maximum peak to peak value and the incident pump power 
respectively. 



in blue, and the numerical integration results are seen in 
dashed red. The two subplots include data obtained with 
E15 (panel («)) and E16 (panel (ii)), respectively. The 
numerical results were calculated using the correspond- 
ing parameters for each device as discussed above. The 
comparison shows a good match between the model and 
the experimental data for both cases. 



VI. ANALYTICAL RESULTS 

Fig. [9] shows a comparison between the measured SM 
frequency and the one predicted by Eq. (P^. for data 
taken with E15 (panel {{)) and E16 (panel (ii)) devices. 
Eq. ([22|) is expected to hold when the input power is 
slightly above the first SM power threshold. Indeed, 
when using an experimentally measured damping rate 
of 7s = 26 MHz for Elfti^ and a fitted damping rate of 
7s = 23.1MHz for E15, the model yields a good agree- 
ment for both devices. 

The spectral power density of the reflected power from 
the resonator is predicted by Eq. ((26l) for the case where 
the pump power is slightly above the first SM power 
threshold, namely regular SM with a rather low frequency 
occurs. This prediction is compared with typical experi- 
mental results, obtained with E15, in the upper subplot 
of Fig. [TOl The noise has a negligible influence on the 
SM characteristics in that power range. On the other 
hand, the dynamics of the system is governed by the noise 
on the edge of the SM, as is demonstrated in Fig[7ja). 
The frequency domain of the numerical results at that 
region is compared with typical experimental results, ob- 
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FIG. 9: (Color Online) SM frequency as a function of the 
normalized injected power i9, measured with E15 (upper) and 
E16 (lower) devices. 
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FIG. 10: (Color Online) The solid blue curves in the upper 
and lower subplots show typical SM experimental results in 
the frequency domain, obtained with E15 and E16 respec- 
tively. The dashed red curve presents the model prediction 
according to Eq. (|26[) . The solid red line was obtained by nu- 
merically integrating the model's equations of motion at the 
first threshold power with nonvanishing noise (Fig[7ja)) and 
evaluating the spectral density. 



tained with E16 in the lower subplot. Both comparisons 
show a good agreement. The numerical results indicate 
that the model predicts a strong noise rise near the SM 
power threshold. Theoretically, such noise amplification 
is expected to increase when the system is approaching 
a threshold of instability, where a linear theory predicts 
an unbounded increase of fluctuation's.^ which only sat- 
urates due to the high order nonlinear terms^. As pre- 
dicted theoretically^ the same mechanism generates large 
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signal amplification, as was indeed observed in RefiM. 

VII. CONCLUSIONS 

We report on a novel nonlinear behavior, where SM 
is generated in SC microwave stripline resonators. This 
phenomenon is robust and occurs with all of our devices, 
despite differences in geometry, and at various resonance 
frequencies in each device. A theoretical model accord- 
ing to which the SM originates by a thermal instability is 
proposed, to account for our findings. In spite of its sim- 
plicity the model exhibits a good quantitative agreement 



with the experimental results. These devices can serve 
as ultra-low noise amplifiers with possible applications in 
the field of quantum data processing. 
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